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Abstract 

The statistics of charge transport across a tunnel junction with energy-dependent scatter- 
ing is investigated. A model with quadratic dispersion relation is discussed in general and, 
independently, in the two limiting cases of a large detector and of a linear dispersion relation. 
The measurement of charge takes place according to various protocols. It is found that, as a 
rule, the statistics is expressed by means of time-ordered current correlators, differing how- 
ever from the conventional ones in the ordering prescription. Nevertheless binomial statistics 
is confirmed in all cases. 

1 Introduction 

The statistics of charge transport AQ through a junction is encoded in its generating function 

X ( A ) = (e*">y = J] ^ <(AQ)»> . (1) 
n=o n - 

Different notions of charge transport, reflecting different measurement protocols, have been pro- 
posed. Each of them has been more or less closely associated to current correlators, which involve 
some time ordering prescription (like the Dyson [l] or Keldysh [2j|5] variants) or do not [6]. The 
resulting statistics differ, but are always expressed in terms of the transparency T of the junction. 
Somewhat unprecisely, the 3rd cumulant is 

«(AQ) 3 »cx - 2T 2 (I - T) , «(AQ) 3 »(XT(I - T)(l - 2T) , (2) 

depending on whether time ordering is forgone, resp. imposed. Matters are however more subtle, 
since it was shown [4j|7] that, in some model with linear dispersion and depending on further 
details, the first result (p| arises even if the correlators are taken to be time ordered. The issue 
was clarified in 18], where it was shown that: (i) as a rule, both variants of time ordering require 
amendment, giving way to unconventional prescriptions; and (ii) in the particular case of that 
same model, the new procedure yields the second result Q, at least in the Dyson variant. 

The purpose of the present work is again twofold: First, to confirm item (i), though from a more 
general perspective, emphasizing the role of gauge coupling; and, second, to compute cumulants on 
the basis of the new procedure for models or variants not considered before. In particular, we will 
consider a model with quadratic dispersion and the cases in which the time ordering prescription 
reduces to the conventional one. In contrast to previous studies, binomial statistics is confirmed 
in all regimes. 
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For concreteness we present our results first within the setting of [9j: A system which is endowed 
with charge and coupled to a spin |, the purpose of which is to serve as a detector, specifically as 
a galvanometer. Let the Hamiltonian of the combined system be 



V H{-\/2) 

where the coupling A and the Hamiltonian H(X) will be specified later. For now we note that the 
spin precesses about the 3-axis, since cr 3 is conserved. 

Let P® pi be the initial joint state of the system and of the spin, and set {A) = tr(PA), where 
A is any operator of the system proper. Then 

X (A) = (jH{-\mt e -i H {\ii)t^ (4) 

is the influence functional describing the spin state alone at a later time t, 

p++ p+ -)^ Pi ={ p+ t x] p+-xW\ (5) 

P-+ P—J \P- + X(-A) P— J 

the representation being again in the eigenbasis of (73. With a grain of salt it may also be identified 
with the generating function in Eq. (|T|), see item C3 below for details. 
We restate Q as 

X(A) = (Vexp (ij di'ffiC-l.o) Texp (-ij dt' H^, if)\ \ , (6) 

where H is the Hamiltonian of the isolated system and 

fli(A,t) = e lHt {H{\) - H)e- lHt 

is the Hamiltonian in the interaction picture; moreover T , T denote the usual and the reversed 
time ordering. We recall that the time ordering is supposed to occur inside the integrals once the 
exponentials are expanded in powers of A. Eq. ^ follows by inserting 1 = exp(— iHt) exp(iHt) in 
the middle of the expectation Q and by performing a Dyson expansion. 

Let Q be the charge to the right of the junction. It is considered to be a primitive quantity, 
corresponding to Q®1 for the combined system, but independent of A. By contrast and as implied 
by charge conservation, the current is then a derived quantity, namely the rate of charge: 

/(Aa 3 /2) = 1 = i [H(Xa 3 /2),Q (g> 1] = ( /(A Q /2) /( _° A/2) ) , (7) 

where /(A) = i[H(X),Q]. 

We shall next present two general types of Hamiltonians H(\) to be used in Eq. (|3j. They 
are constructed from the Hamiltonian H and from the charge Q. At first our goal is to emphasize 
structure. Later a more concrete physical meaning will be attached to the construction by means 
of a series of examples for H and of remarks about H(Xas/2). 

The general Hamiltonians H (A) are as follows. 

Al. Linear coupling. The Hamiltonian is 

H{\) = H-XI, 

where / = i[H, Q] is the bare current through the junction, i.e. in absence of coupling. 
Hence, 

Hj(X,t) = -XI(t) , 
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where I(t) is the bare current in the interaction picture, and Eq. (|6]) reads 
X(A) = / a* exp fif J dt' I(t')\ T exp (i\ J dt' I{t') 



(8) 



which is the announced and well-known relation with current correlators (Keldysh time or- 
dering). This setting is closely related to that of |10] . 

A2. Gauge coupling. Gauge transformations are generated by local charges Q. The Hamiltonian 
H transforms as 

H(X) = e lXQ He- lXQ (9) 
= H-XI-i^-[Q,I] + 0(X 3 ), (A-0), (10) 

where I = i [H, Q] . (A specific model illustrating that coupling scheme on a spin will be 
considered shortly in CI.) Local currents are obtained by varying the gauge, 

H'(X) = ±H(\) = i [H(X),Q] = -/(A) , (11) 

which results in 

/(A) = e lXQ Ie~ lXQ . 

The current parallels the kinematic momentum of a particle in that its value is gauge invari- 
ant, but its representation as an operator is not, i.e. /(A) ^ J as a rule. This is reflected in 
the current correlators, which are now those of 



ffi(A,t) = - f d\'I(X',t), 
Jo 



showing that Eq. (|8| is in need of amendment, and not just by replacing J(t') with /(+A/2, t'). 
The central result is that x(A) can still be expressed by bare current correlators; however 
Eq. ^ has to be modified as 

x(A) = (V^xp(zf |V/(0) F'expfif J" (#1(0)) , (12) 

1 — * 

where the unconventional ordering T means that the derivative in 

dt 

has to be taken after the time ordering, 

T* (/(t x ) • • • 7(t„)) :=±...±T (Q(ti) ■ ■ ■ Q(tn)) , (13) 
(Matthews' time ordering [ll]); likewise for 2* . Eq. (12 1 follows from ^ by the identity 



f exp f-i J dt' Hi(X,t')^j = ^T* exp (iX J dt' J(t') 



(14) 



which will be discussed in Sect. [5j In summary: The T-ordering in connection with Hi is 
equivalent to T* -ordering in connection with —A J, see Eqs. (|6 12 1. 
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Item Al is included in A2 as the special case in which 



[Q,J] = 0, (15) 

see Eq. (10 1. Then the T and T*-orderings of currents may be used interchangeably. Moreover, 
the operator of current ^ becomes insensitive to the inclusion of spin: /(A) = /. 

Fairly concrete examples illustrating the general types are provided by a single particle moving 
on the line. Depending on its position x e R, the charge Q to the right of the junction is 1 or 0, 
and is in fact implemented as the multiplication operator by the Heaviside function 6(x). That 
function can be replaced quite conveniently by a smooth version thereof, Q = Q(x); the size of 
the region where they differ may be loosely associated with that of the detector. Through second 
quantization the examples implicitly describe many independent particles, too. 

Bl. Linear dispersion. The Hamiltonian H = p+V(x) on L 2 (M) describes a right moving particle. 



The current / = i[H, Q] = Q'{x) satisfies \Q,I\ = 0, i.e. (15 1, whence the T-ordering suffices 
as a rule. By combining two copies of such a model, scattering between channels of left and 
right movers can be included. However Eq. ( |15[ ) fails in the limiting but computationally 
simple case of Q(x) = 9(x) and V(x) a point interaction at x = [81. As a result, T*-ordering 
is required in this special case. 

B2. Quadratic dispersion. The Hamiltonian is H = p 2 + V(x). Then 

I = i[H,Q] = P Q'(x) + Q'(x)p, (16) 
[Q,I] = 2iQ\x) 2 # 0, (17) 

which calls for T*-ordering as a rule. However, in the limiting case of an ever smoother 
transition function we have Q' 2 — * (in any reasonable norm). As a result, T-ordering 
should suffice in that special case, as we will confirm. 

A few remarks will now address the role of the spin. 

CI. We recall [12] that the Hamiltonian (|3] [oj is physically realized by a spin coupled to the 
current flowing in a wire, as we presently explain. Let the straight wire run along the 1-axis 
and let xq be the position of the spin in the 12-plane. The vector potential due to the spin 
a/2 is 

where f(x) = n\x — a?o | — 1 ■ (More general functions / are obtained by smearing the position 
of the spin.) A particle in the wire couples to the spin through (p — A) ■ e\ = p — A - ef, where 

A ■ ei = (ex a V/) • - = (e x a V/) • e 3 — = — — , 
2 2 0x2 2 

since by the stated geometry e\ a V/ lies in the 3-direction. We note that df/dx2 = 0(r~ 2 ), 
(r = \x\ — > oo); hence, along the wire, df/dx2 = XQ'(x) for some coupling A and a function 
Q(x) of the kind described above. In summary: As p gets replaced by 

(p- A) ■ S x = p- XQ'(x)^ = e ^Q^ pe -^Q^ ; 

so does H by H(Xa 3 /2). 

C2. In the previous item AQ'(x) arises as a connection. It appears with the replacement A — » 
A173/2, by which it acts non-trivially on the spin. As we presently explain, that property pro- 
vides a geometric mechanism (though different from the physical one) by which the rotation 
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of the spin becomes a counter of the transported charge. The mechanism somehow resembles 
that of a screw, whose motion rigidly links rotation with translation. More precisely, a state 
ip of the combined system obeys parallel transport along the line if 



or equivalently if 



(p-AQ'(aOy)^ = 0, 



— = iXQ (x)— - tp. 
ax 2 



Given that a spin state ip changes by dtp = ~ i(a- e/2)ipd0 under a rotation by dd about e, the 
condition states that a charge transport dQ = Q'(x)dx is linked to a precession dQ = — AdQ. 

C3. The conclusion of the previous item can not be reached for the physical evolution of the 
combined system, as generated by the Hamiltonian ^ , at least not without further ado. In 
Bohr's spirit |13| and in elaboration of |12] it is convenient to talk about the apparatus in 
classical terms, here a classical spin s e K a , (\s\ = 1). The spin components Sj have Poisson 
brackets {si,Sj} = tijkSk and the Hamiltonian is H{\s^). In particular Eq. ^ would be 
recovered by quantization. The equations of motion are 

Si = {si,H(Xs 3 )} = Ai/'(As 3 ){si,S3} , 

or, by (fill 



s = — A/(As3)e*3 a s. 

The angle of precession thus is 

6 = -Xq, (18) 

revealing the charge q that has flowed during a time interval [0, t]. In this context \ can be 
interpreted as the generating function of the transported charge: 



X (A)= \dqx(q)e iXq . 



Indeed, since 

n (ff\- ( P++ P+- er 
is the state p\ of the spin after precession by the angle 6, its final state §5§ is (9] 



I 



Pi = dqx(q)Ps(-M) 



In view of ( |18[ ), this is consistent with x(l) being the probability (density) of transport q. 
as claimed. The interpretation is however hampered by the fact that \(q) may fail to be 
positive 1 14 . 



We should also mention an earlier approach (T 15 to charge transport, which does not explicitly 



model a detector. It is based on two measurements of the charge Q, occurring at times and t. 
The transported charge is then identified with the difference AQ of their outcomes. The associated 
generating function is 

X(X) = (e^Q(t) e - lX ^ , (19) 

at least in the case when the initial state p is an eigenstate of Q or an incoherent superposition of 
such, i.e., for [p, Q] = 0; then x actually agrees with the expression Q. We will use the definition 
beyond this restriction, because it is irrelevant in the limit of large times. (See however 1 16] for 
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the unrestricted definition.) By (e tHt e lX( ^e %m )e lX< ^ = e im (e lXC ^e lHt e lXQ ) we may restate the 
generating function in a form closer to Q . 



e iHt e -iH(\)t 



with H{\) as in Eq. ([9|; and further in terms of current correlators by means of (14 1 

*(A) = : exp(i\j I(t') 



(20) 



where the star can again be dropped under the assumption (15 1. 

We shall now describe the main result. It confirms the binomial statistics of charge transport 
in a variety of situations. Specifically, in the long-time limit the 2nd and 3rd cumulants of charge 
transport are 



ton k(AQ) 2 » = ~ p dE T(E) (1 - T{E)) , 

ton k(AQ) 3 » = i ^ dE T{E) (1 - T(E)) (1 - 2 T(E)) , 

•'Mr. 



(21) 
(22) 



where < fi h are the Fermi energies of the states incoming from the right and the left sides of 
the junction, and T(E) is the transmission probability (transparency) at energy E. (Eq. pTj ) was 
first obtained in [IT] without any time ordering prescription.) The results apply to 



Dl. either generating function, Eq. (12) or ( 20 I ; 

D2. Hamiltonians with linear or quadratic dispersion relation (see items B, but in second quan- 
tization) ; 

D3. independently of how sharp the jump of Q is, i.e. of the width over which Q(x) differs from 
0(x). 



Of some interest is the way that independence arises. The time ordering (13 1 is spelled out for 
n = 2 as 



T* (1(h) i(t 2 )) = (Q(h)Q(h)9(h - h) + Q(h)Q(h)0(h - h)) 

oh oh 

= T(i(h)i(h)) + [Q(h),i(h)]S(h - h) , 



(23) 



and thus differs from the usual one, T, by contact terms supported at coinciding times; likewise 
for n = 3, where 

^r*(hhh) = ^T(hhh) + 3<5(i 2 - hfT(h[Q 2 ,h]) + S(h - h)S(h - t 3 )[Qi, [Qi,h]] (24) 

with the shorthand notation li = 1(h) (for general n, see [8]). Depending on circumstances, the 
terms in the expansions contribute variably to the invariable results (21 22 1, as we now detail. 

Items D2, D3 come with interpolating parameters: The Fermi wavelength Af, with Ap 
the linear dispersion is approached through rescalings 



as 



±P 



of right and left movers; and the width I of the transition region. In the limit I — * go [18] of an ever 
larger detector Q ceases to be defined. In the opposite limit I — »■ and in the case of quadratic 
dispersion (Af > 0), the commutator [Q,I] diverges even in the sense of distributions, because of 



(> 



A F = oo A - 



A F = o 



T 



1 = T 



I = oo 



Figure 1: Matthews' vs. usual time ordering of current correlators. The parameter range (i,Ap) 
of models is shown as a square, which includes solid parts of the boundary. In the limits (thick 
arrows) of linear dispersion (Xp — » 0), or large detectors (I — * oo), the contact terms appearing in 
Matthews' time ordering, T*. vanish. In these limiting cases there is agreement with usual time 
ordering, T. 



Q'(x) —* 6(x) (see item B2). One can though discuss the limits of the correlators (21 22) in these 
limits, and of their parts, but not those of the models themselves. However a model with Af 
I = exists [8], describing a scatterer and a detector which are both pointlike and coincident 



0. 



Let us discuss Eq. (|20| first. The contributions of contact terms to the cumulants (21 22 1 are 

as shown in Fig. [I] As for Eq. ( 12 1 the same 



oo, 



non-trivial, except in the limits for which I/Xf 
is true for the 3rd cumulant: however for the 2nd the contact terms cancel between T and T 



Relation with the literature is eased through books and the review articles on noise and count- 
ing statistics, and among them (7|[l4}[l9j. 



The plan of the article is as follows. In Sect. [2] we introduce a model with quadratic dispersion 
relation and review the main features of a limiting case with linear dispersion. In Sect. [3] we will 
explain the broad structure of the computation of the cumulants and emphasize the methods. The 
first part of Sect. [4] is devoted to the detailed derivation of asymptotic binomial statistics for the 
model with quadratic dispersion. In the second and third parts, the limiting cases of a large detec- 
tor and of linear dispersion are given independent derivations. In Sect. [5] we recall the reason for 
the T* time ordering and discuss the equivalence between the methods used here, resp. elsewhere, 
as e.g. in |4|. Finally, two appendices collect some auxiliary results. Appendix [A| contains the 
long-time limits of some distributions, while in Appendix [B] matrix elements of charge and current 
are computed. 



2 The models 

2.1 The quadratic dispersion model 

We consider two conducting leads connected through a junction and model the whole device by 
independent fermions moving on the real line. The single-particle Hamiltonian is H = p 2 + V(x), 
acting on the Hilbert space L 2 (JBL). The kinetic energy is quadratic in the momentum p = —id/dx; 
the potential V describes the junction and vanishes away from it, i.e., outside of some interval 
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[— xq, xq\. The potential will enter the discussion only through its reflection and transmission 
amplitudes, r(k) and t(k). They can be read off from the Lippmann-Schwinger (LS) states \ij)k)'. 
Continuum eigenstates of H of incoming momentum k and eigenvalue E = k 2 have wave- 
functions ipfc given outside that interval as 



k > : ip k (x) 
k < : t/j k (x) 



e lkx + r (J.)e lkx , (x < -X ) 

t(k)e tkx , (x > x ) 

t(k)e tkx , < -x ) 



Note that states with k > (k < 0) have incoming parts that are right (left) moving. By the 
Schrodinger equation the scattering matrix 

~ \r(k) t(-k) 

is unitary. In particular, the transmission and reflection probabilities are even in k, whence 
T(fc 2 ) := \t(±k)\ 2 and R(k 2 ) := \r{±k)\ 2 , and satisfy 

T(E) + R(E) = 1 . (26) 

LS states form a (continuum) basis of L 2 (M.) normalized as 

ifd* 1^X^1 = 1- (27) 

^ JR 

Time-reversal invariance of H is, incidentally, a property which is not relied upon, in that the 
above discussion still applies when p is replaced by p — A(x), at least as long as A has the same 
support properties as V. 

As mentioned in the introduction, the charge to the right of the junction may be implemented 
on L 2 (M.) as a multiplication operator, Q = Q(x). More specifically, we assume 

{ ; <xo \ (28) 

^1, (X » Xq) . 

The left and right leads are assumed to be reservoirs with energy levels occupied up to Fermi 
energies /i L , [i R > biased by V = fi L — /U R > 0. The occupation of LS states thus is 

p(*) = | 1 ' ( 29 ) 
I , otherwise 

where fc LjR = (/x LiR ) 1,/2 ; or for short p(k) = 6{k e J) where J := [— fc R , fcj. More precisely, p = p(k) 
is the single-particle density matrix (0 < p = p* < 1) of the many-particle state (•); actually p 
determines a quasi-free fermionic state (•), which for practical purposes means that expectations 
of many-particle operators can be computed by means of Wick's rule. As a matter of fact, for p a 
projection as in (29 1 the state (•) is necessarily quasi-free. 



2.2 The linear dispersion model 

We briefly review the main features of the linear dispersion relation model used in [8] and underlying 
the computations of Sect. |4.3| For a more detailed exposition, we refer to the original paper. 

In the limit of long times (or low frequencies) it appears appropriate to linearize the dispersion 
relation near the Fermi energy. A suitable model arises by reinterpreting the two leads on either 
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side of the junction: Rather than viewing them as non-chiral half-lines, they are now (full) chiral 
lines. In absence of scattering, which now amounts to a cut junction, the Hamiltonian is linear in 
the momentum p = —id/dx and is given as 



Ho 



p 
p 



on i 2 (K) © L 2 (IR). A point scatterer is then placed at x = 0; it results in a unitary scattering 
matrix 



s = 

which is independendent of energy. In particular, T = \t\ 2 = \t'\ 2 and R = \r\ 2 = \r'\ 2 still satisfy 



Eq. (26 1. The single-particle charge operator is the projection onto the right lead, 




1 



and the initial single-particle density matrix is the projection 

P 

representing two infinitely deep Fermi seas biased by V = fi L — p R > 0. The condition [p, Q] = 0, 



(O(^-p) 

V o 0(mk-p) 



underlying the unrestricted use of the generating function (19 1, is satisfied here. 

A feature of the model is that the scattering process is instantaneous in the sense that the 
position of the point scatterer coincides with that of the detector. As a result, [Q, I] # ( (8], 
Eq. (3.17)), and the contact terms arising from T*-ordering matter. In terms of the discussion 
given at the end of the introduction, the model has vanishing length scale I = 0. 

However the scattering process can be regarded as strictly causal by separating the two positions 
by I > 0. This is achieved by replacing the charge operator Q by its regularization Qi := Q0(|x| > 
I), and accordingly the current / by // := i[H, Qi] = Q[S(x — I) — 8(x + I)]. Then the commutator 
[Ql, I{] vanishes and with it all the contact terms. 



3 Overview 



Before engaging in the detailed computation of the cumulants (21 22) it is worthwhile giving an 
overview of the methods involved, and illustrating them in simple instances. The physical setting 
has been discussed at length in the introduction and will be recalled only briefly. We consider two 
leads separated by a tunnel junction, with particles in an initial multi-particle state (•). We investi- 
gate the statistics of charge transport, AQ, across the junction and during a time t. Specifically, we 
are interested in its moments ((AQ)™), determined as the expansion coefficients of some generat- 
ing function, see Eq. Q; and actually in the long-time limit of the associated cumulants (((AQ)")). 



Generating functions. In the introduction two distinct generating functions were presented: 



(30) 



X(A) = (V* exp dt' I{t') 

X(A) = (V* exp ^| J* dt' I(t')^j T* exp ^| J* dt' I{t') 
where 

I(t) = e lHt Ie- lHt (31) 
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is the current across the junction. It is expressed in terms of the charge Q(t) to its right as 
I(t) = dQ(t)/dt, whence 

I = i[H,Q]. (32) 
We shall refer to \ an d x as the generating functions of the first and of the second kind, respectively. 

Results. The 2nd and 3rd cumulants exhibit asymptotic binomial behavior. 

ton k(AQ) 2 » = ~ p dE T{E) (1 - T{E)) , (33) 
ton k(AQ) 3 »= i \^ dET{E){l-T{E)){l-2T{E)) , (34) 

*-»CO t 27T l„ 

where T(i£) is the transparency and /j, LjR are the Fermi energies on the left and right leads, in 
various instances and for either generating function. Specifically: 

- in the quadratic dispersion model, contact terms matter up to the special case of an ever 
smoother step of the charge operator Q(x) (see Section 2.1|. 

- in the linear dispersion model, contact terms vanish, except in the special case of instanta- 
neous scattering (see Section 2.2 1. 

In the rest of this section we address the methods used to obtain the results from the generating 
functions. 

< — * 

T*-ordering. It is convenient to recall the expansion in contact terms for T -ordered products. 
With the shorthand notation A { = Afa), A = I,Q, Eqs. (|23]|24]) read 

^*{hh) = 1 T{hh) + S{ti-t 2 )[Quh], (35) 
T\hhh) = Tihhh) + 3S(t 2 - t^Tih^h]) + 6(1, - t 2 )S(t 2 - i 3 )[Q x , [Qx,h]] . (36) 

The ordering by 2* yields the same expansions, up to a minus sign for contact terms involving 
an odd number of commutators. A general expression for products of all orders may be found in [8] . 

Cumulants. Based on the generating function \ °f the first kind we have 

«(AQ) 2 »= fd 2 t«2-*(/ 1 / 2 )»= fd 2 i«7W 2 )»+ fdtx «[<?!,/!]», (37) 
Jo Jo Jo 

«(AQ) 3 »= f d 3 i«T*(W 3 )» 
Jo 

= f &H0{hI 2 h))) + z[ d 2 t«^(/ 1 [Q 2 ,/ 2 ])»+ f dii«[Qi,[Oi,/i]]», (38) 
Jo Jo Jo 

where d n t = dt\ ...dt n . At first, similar equations are obtained for the moments by means of 
Eqs. (30 1 and (35 36). Moments can then be replaced by cumulants; indeed, their combinatorial 
relation is universal, and hence the same on both sides of the equations. 
Based on the generating function x of the second kind we likewise find 

«(AQ) 2 »= i£d 2 n«T*(/ 1 / 2 )» + 2«/ 1 / 2 » + «r*(/ 1 / 2 )»] = £d 2 t«/!/ 2 » (39) 

«(AQ) 3 » = i £ d 3 i [«7**(W 3 )» + 3«2**(IiJa)l3» + Hh^*(hh))) + «^*(W 3 )»] , 

(40) 
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where the 3rd cumulant may also be expanded using (35 36) for both time arrows. We note 
that the cumulants of the second kind involve -ordered current correlators already present in 
those of the first kind, and more. However, due to the symmetry between usual and reversed time 



orderings, the contact terms in the 2nd cumulant (39) mutually cancel 



Wick's rule. The many-particle state (•) is the quasi-free state 20 determined by the single- 
particle density matrix p. Let A be the second quantization of the single-particle operator A. 
Correlators of second quantized operators can be reduced to the level of first quantization thanks 
to Wick's rule. In particular, with p' := 1 — p we have 

((A)) = (A) = 0, (41) 

«A§» = tr(pAp'B) , (42) 

((ABC)) = tr(pAp'Bp'C) - tx(pAp'CpB) . (43) 

The expressions follow from the usual formulation of the rule which involves creation and an- 
nihilation operators %j)*{a), ip(b) of single-particle states a, b: Expectations of products of such 
are computed by way of complete contraction schemes and reduced to just two kinds of contrac- 
tions, (ip* (a)tp(b)) = (b\p\a), (i/;(b)ip* (a)) = (b\p'\a), with further ones vanishing. The second 
quantization A i— > A is defined for rank-one operators A = \ai)(a,2\ as 

A = i/j*(a 1 )ip(a 2 ) - (a x \p\a 2 ) 
and then extended by linearity in A. We stress the "zero-point subtraction" of (ai\p\a-z) = tr(pA), 



which implies (A) = 0, but drops out from higher cumulants. The l.h.s. of Eq. (42 1 gives rise to 
just one non-vanishing connected contraction scheme and thus equals 

«^*( ai )^(a 2 )V*(6i)^(fo 2 )» = <^>iM& 2 )X^(aa)V'*(&i)> = (MpKX^Ip'I&i) , 

in agreement with the r.h.s.. 

The charge and current operators mentioned earlier in this section are meant in second quan- 
tization. We will henceforth denote them by Q and /. 

GNS space and Schwinger terms. The reader may skip this item. It is in fact about some 
fine points which remain without practical consequences. Strictly, ip*(a), ip(b) act on the GNS 
space |20| of (•) but, as we explain below, one may sometimes pretend it is replaced by Fock space. 
We only consider the case when p is a projection. An operator A admits a second quantization A 



if B = [p,A] is Hilbert-Schmidt, i.e. tr(B*B) < go. The traces (42 43 1 exist if A and the other 
observables satisfy that condition. The (first quantized) operators Q and I do so in both models 
of Sect.[2j However, ti(pl) is well-defined only in the model with quadratic dispersion, and tr(pQ) 
in neither. 

For p = the stated condition becomes trivial, and the GNS space is the Fock space. If its 



operators A are used on another quasi-free state p, then (42 43 1 are still valid if the traces exist 



but (41 ) is to be replaced by ((A)) = (A) = tr(pA). The difference is in the "zero-point subtraction" 
which may diverge, even for [p, A] Hilbert-Schmidt. The point of the GNS space is that A still 
remains defined there. 

On the GNS space we have [20] 

[A, B] = [AJ3] + S(A, B) ■ 1 , S(A, B) = ti(pAp'B) - tr(pBp'A) , 
where the last term is known as a Schwinger term. We infer 

A{t) =Mt) + i[dt' S{H, A(t'))l , (44) 
Jo 

(([A,B])) = S(A,B). (45) 
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Upon pretending that the GNS space is just Fock space, we have [A, B] = [A, B] and A(t) = A(t). 
However, Eq. (451 still holds true, because (42) does. The same conclusion is obtained from 



(([A,B])) = (([A,B])) = tr(p[A,B]). 

Let us comment on the significance of Schwinger terms for the cumulants (37 38 1 , where /j 
is now to be read as I(tA (and likewise for Q.A. First, it may be replaced by I(ti), because the 
difference seen in (44 1 drops out from the results. Second, the contact terms <<([Q(£i), I(£i)])) and 
(([<3(ii), [Q(ti), /(£i)]])) are properly Schwinger terms. Informally however they may be under- 
stood in the context of Fock space, as discussed. 



A simple case. We illustrate the methods by considering a simple example: The 2nd cumulant 



of the second kind, Eq. (39 1, for the model with quadratic dispersion and with state (29 1. Traces 



may be evaluated using the basis (27 1 of LS states. We so obtain 



«li/ 2 » = tr{php'h) = jd 2 k(l\ph\2X2\p'I 2 \X) = 7^)2 J d 2 fc e ^-^H^)|<l|/|2>| 2 , 

R 2 JxB\J 

(46) 

with the shorthand notations | i) =| V'fc;), Ei = kf, and J = [— k R ,k L ]- The last equality follows 
from Eq. (31 I and the eigenvalue equation H \ i) = Ei \ i). 



Time integrals. The long-time limit of (|39| now calls for 
im - f 



hm - J* d 2 t e^~ E ^-^ = 2nS(E 1 - E 2 ) 



2tt 
2;*i 



{8{kx-k 2 )+5{k l + k 2 )). 



The first equality is by Eq. 



below and the second by Ei 



Only one of the diagonals 



k\ = +k 2 openly intersects the integration domain J x R\J in Eq. (46), and in fact just for 
ki = —k, 2 e [fc R ,fc L ], see Fig. [2} Hence 



\ r 






x -1 



Figure 2: Integration over k\, k 2 in Eq. {J^6). The integration domain J x R\J (shaded) and its 
intersection with E\ = E 2 (diagonals). 



hm J«(AQ) 2 » 

i->oo t 



\ ki I<1|/|-1>| 2 



where |— 1) :=\ip_k i y. Appendix [A] collects further long-time limits of distributions. 



(47) 
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Matrix elements of current. The current / is given in Eq. (16 1. We compute the relevant 
matrix element for k\ > and observe that Eq. (25 1 distinguishes between the cases +k > 0; 
however only the expressions for x > xq matter here because of the support properties of Q' seen 
in (j28j: 

<1|/|-1> = <l|pQ'(;r) + Q'(x)p\-1) 

rOD 

dxt(hje~ tklX ( P Q'(x) + Q'(x)p)(e-' lklX + r(-fc 1 )e I ' £lX ) 

J— (X) 

= 2k 1 t(k 1 )r(-k 1 ) dx Q'(x) = 2k 1 t(k 1 )r(-k 1 ) , 

J— CO 

where the third equality is by partial integration. Hence 

|<1|/|-1>| 2 = {2k 1 YT{k\)R{kl) = {2k l fT{k\){l T{k\)) , 



(48) 



(49) 



and Eq. (33 1 follows by substituting E = k\ in Eq. (47 1. In the simple case considered we thus 
confirm the binomial statistics. 

Further computations of matrix elements of current may be found in Appendix[B] In particular, 
we mention 

<1|I|1> = 2k 1 T(k 1 ) , (ki > 0) . (50) 

Matrix elements of charge. The cumulant (39 1 we just computed is the simplest among ( 37|40 1 
in that it does not involve the charge operator Q. In preparation of the other cases, it pays to look 
at the relation between matrix elements of Q and of /, still though within the model of quadratic 
dispersion. In view of the support property (28) of Q(x), computing its Fourier transform demands 
a regularization at x — > +oo: 



rOD 

Q(k) = lim dx Q(x)e-^ k -' l ^ x = (-») lim 

e|0 eiO 



Q'(k) _ Q'(k) 
k-is [ l 'k-iO 



(51) 



The result is a distribution in k, and so is (1|Q|2) in fc 2 ; whereas (1|/|2) is a smooth function of 
these variables. By (32 1 we have the equation 



<l|/|2> = i( J B 1 - J B 2 Xl|g|2>, 



(52) 



which however can not be uniquely solved for (1|Q|2), because the distributional equation xF(x) = 
admits the non-trivial solutions F(x)cc5(x). In fact, in view of the Sokhatsky-Weierstrass (SW) 
formula 

' ' 2mS(x), (53) 



the general solution is 



i<l|Q|2> 



x — iO x + iO 
<1|7|2>(+) 



+ 



<1|7|2>(- 



(54) 



E 1 -E 2 + i0 E x - E 2 - i 

where <1|/|2> = <1|/|2>W + <1|J|2>^) is any split of the matrix element of current. It takes <1|Q|2> 
to make it unique, at least up to terms vanishing for Ei = E 2 , which may still be shifted between 
the two terms <l|Zj2>t ± >. 

In summary: An expression for (1 1 1 \ 2) does not entail one for (1 1 Q \ 2); rather conversely, 
including the split. Such expressions will be derived in Appendix [B] An impo rtan t case is when 
k\ ^ k 2 , whence E\ = E 2 arises by k\ = —k 2 \ then (1|/|2)( - ) = (1|/|2) and Eq. ( 54 1 simply reads 



Another case is 



<1|Q|2> = (-i) 
<-l|/|-l>(-) 



<1W2> 



E x - E 2 - i 
2k x R{k x ) , (fci > 0) . 



(55) 
(56) 
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4 Derivations 



4.1 The quadratic dispersion model 

Using the methods introduced in the previous section we shall derive the asymptotic binomial 
distribution (33 34) for both generating functions, \ an d X- We shall do so first for the model 



with quadratic dispersion relation of Section |2.1| It will become evident that contact terms are 
crucial. In other words Matthews' time-ordering can not be replaced by ordinary time-ordering, 
except in limiting cases, if the correct result is to be found. Two such cases, namely that of a large 
detector and of a linear dispersion, will be given independent treatments in the following sections. 



2nd cumulant of the first kind. The cumulant is given in Eq. (37 1 as 

«(AQ) 2 »= f d 2 t«T*(/ 1 f 2 )» = A + B 



with 



A 
B 



f d 2 t«^(f 1 f 2 )» 
Jo 

f dti«[Qi,/i]» 
Jo 



! fdtif 
Jo Jo 



dt 2 «/ x / 2 » 



(main term) , 
(contact term) . 



The connected correlators are computed by Wick's rule (42 1 and the resulting traces evaluated in 
the basis of LS states (25 1. We obtain 

2 



(2tt)2 



i (E 1 -E 2 )(t 1 -t 2 ) 



I<1|^|2>| 2 



B 



{2-Kf 



f d 2 k fdtif df : 
J Jo Jo 

kM\J 

J d 2 M<lM2><2|/|l>-<l|/|2><2|Q|l». 



(57) 
(58) 



Jxl 



In relation to the overview above we shall next (i) discuss time integrals and (ii) express matrix 
elements of charge in terms of those of current. We will do likewise later for all cumulants. In the 
present case the first item concerns only the main term, the second only the contact term. 

(i) The asymptotic long-time behavior of the main term is given by Eq. ( 87 I with x = E\ — E% 
and the substitution t 2 ^> t\ — t 2 : 



dt 2 e i{El ~ E2){tl ~ t2) 



t—+co Ei- E 2 + i0' 



(59) 



as distributions in k\ and k 2 . In Eq. (57 1 (1|J|2) then qualifies as a test function, being essentially 
the Fourier transform of the compactly supported function Q' (x). 



(ii) Within the integration domain (58 1 (1|Q|2) is given by Eq. (55), and (2|Q|1) is then 
obtained by exchanging fej and k 2l or by complex conjugation. 
Collecting terms we so obtain 



2i 



Jx«\J 

dki 

) fe R 



E 1 - E 2 + i 
I<1|^|-1>| 2 



+ 



Ei - E 2 - i E 2 - E x - i 



)I<1|W 



p f L dET(E)(l-T(E)), 
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as claimed. The second equality is by the SW formula (53 1 and the earlier remark restricting ki 
to [fc R , fc L ] (see Fig. pj|; the last one by Eq. (49 1 with E =~kf. 

3rd cumulant of the first kind. Though longer, the computation of the 3rd cumulant retains 
the same two key ingredients: (i) the evaluation of time integrals and (ii) the expression of matrix 
elements of charge in terms of those of current. By Eq. (38 1 we have 

«(AQ) 3 » = f d 3 i <<^*( W 3 )>> = A + B + C 
Jo 

with 

A = f d 3 t ((Tihhh))) = 6 fdixf d* 2 f d* 3 «/i/2/s» (main term) , 

Jo Jo Jo Jo 

B = 3 f d 2 t ((T(h [Q 2 , 7 2 ])» = 3 f dti f dt 2 «£ [Q 2 , 7 2 ] + [Qi, h]h)) (1st contact term) , 
Jo Jo Jo 

C= f dti«[Qi,[Qi,fi]]» (2nd contact term) . 

Jo 

Here Wick's rule d43| is appropriate, whence each term splits into two. 



(a) Main term, (i) We get 

«7i7 2 7 3 » = tr(pI lP 7 2P 7 3 ) - tv(pI lP 'I 3 pI 2 ) 

\^ f d 3 fce 4 ( £l - B2 )* 1 e t(£;2 - £;3)t2 e l(£;3 - El)t3 <l|/|2X2|/|3X3|7|l> 

ivy j 



(2tt)= 

JxR\JxR\J 

- J d 3 /fce l(iJl " iS2)tl e l(E2 - B3) * 3 e l(£3 - £l)t2 <l|/|2><2|/|3X3|/|l>. 

JxR\Jx J 

The exponentials of the second term are obtained from those of the first one by exchanging E 3 — Ei 
by Ei — .E3, which leaves the sum E 2 — E\ unaffected. This symbolic transformation rule may be 
deduced from Wick's rule (43 1. For possibly distinct observables and for the terms as a whole it 
reads: 

- exchange <3|-|1> by <2|-|3>; 

- exchange E 3 — E\ by E 2 — E 3 ; 

- in the integration domain, replace fc 3 e R\J by k 3 e J: 

- apply an overall minus sign. 

In the present case the first item leaves the integrand unchanged. We shall denote the transforma- 
tion by 7(23), as ^ essentially arises by exchanging positions 2 and 3 in the product of operators. 

The time integrals are given by Eq. ( 89 ) with x = E\ — E 2 and y = E 2 — E 3 resp. y = E%— E\. 
It yields 

lim - [ dH ((^(hhlz))) = A 1 + A„ 



with 



<1|/|2><2|7|3><3|/|1> 

JxR\JxI 



Ai = -— ^ d 3 *; 



1 



(2tt) 3 J (E 1 -E 2 + iO)(E 1 -E 3 + iO)' 

<1|7|2><2|7|3X3|/|1> (60) 



An = 7(23)[Ai] = ^3 J" d 3 fc 



(E 1 -E 2 + iO)(E 3 -E 2 + iO) 

Jx«\JxJ 
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(b) Contact terms, (i) The long-time behavior of the 1st contact term is given by Eq. (87 1, 
while the integrand of the 2nd is time-independent. Expanding the commutators we obtain 



lim j f d 2 t«^(f 1 [Q 2 ,i2])» = B I + B„, 

:-+00 t Jo 

lim if dt 1 «[Q 1 ,[Q 1 ,/ 1 ]]» = C I + C„, 

t^+oo t J 



with 



B i = — — J d'L- , 



/ <1|J|2X2|Q|3X3|J|1> - <1|/|2><2|/|3><3|Q|1> 
I_ (27r)3 J \ E 1 -E 2 +iO 



JxR\JxR\J 

<1|Q|2)<2|/|3)<3|I|1) - <l|/|2)<2|Q|3)<3|I|m 
Ex~E 3 +iO J ' 

Ci = J d 3 fc «1|0|2><2|Q|3><3|/|1> - 2<1|Q|2><2|/|3><3|Q|1> + <1|/|2><2|Q|3><3|Q|1» . 

JxR\JxM\J 

and B n = 7^3) Pi], C n = 7^ 23 )[Ci]. 

(ii) We distinguish between matrix elements (i\Q\j} as to whether ki and kj belong to the same 
or to different sets among J and R\J. In the first instance Eq. (55) applies, whereas in the second 



its generalization (|54j) is required. For example the first two terms in the integrand of Bi become 
<1|/|2><2|/|3>( S )<3|/|1> <1|/|2><2|/|3>«<3|/|1> 



s 



(E 1 - E 2 + i0) (E 2 -E 3 + si0) {E x - E 2 + iO)(E 3 - E l - i 0) 



where in the second term we used the identity <2 1 1 | 3> = (2\I\ 3>W + <2 1 1 \ 3>H . In view of 
the integration domains the splitting will more generally affects (2|7|3) in a\ and (3|/|1) in an, 
(a = A, B, C). As a result each term aiji may be written as 



°i=(2^8 \ d3fc E a i (S) (^- J Bi^2-^)<l|/|2><2|/|3>( s )<3|/|l>, 

JxR\JxR\J S = ± 

«ii = ^ \ d3fc S a " (S) ^ 3 - E ^ E * ^)<1|/|2><2|/|3><3|/|1>( S ) 



JxS\JxJ 

for some distributions ai^-K We remark that au^(E 3 —Ex,E 2 —E 3 ) = —av-\E 2 —E 3 ,E 3 —Ex), 
because the matrix elements carrying the superscript (s) in the two cases are also those exchanged 
by the transformation rule 7(23)- Moreover, the dependence of on s is of the form 

(±)rtp tp t? em ai(E 3 - Ei, E 2 - E 3 ) aj(E 3 - E 1 ,E 2 - E 3 ) 
ai {Es El ^ 2 E3) = El -E 3 + i + E 2 -E 3±i ' (61) 

where a\ and ot\ (a = A, B, C) are as follows: 

~ M = - El -E 2 + l > Xl = °' 

Bi = + , Bi = , (62) 

Ex - E 2 + i E 1 -E 2 -i0 E x -E 2 + i0 E x -E 3 + i0 y ' 

ft~ „ I .. . c — _ I 1 



Ex-E 2 -i0 Ex - E 2 - i J5 X - E 3 + i 
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The claim we are heading to is 

Ai + B : + Ci = ~ r dET(E)R(E) 2 , A„ + B„ + C„ = ~ P dE T{E) 2 R{E) , (63) 

which leads to binomial statistics Q in view of TR 2 - T 2 R = T(l - T)(l - 2T). To establish it, 
we observe that the sum, 

A X + B X + Ci = ^3 \ d 3 fc J] A I ( S )<1|/|2><2|/|3>^)<3|/|1> , (64) 

JxR\JxR\J S_± 



likewise involves distributions of the form (61 1 with 

Ai = 2mS(E 1 -Ei)- \— — , Aj = -2m5(Eh. -E 2 ) + - El ~ h " 



E 1 -E 2 + i0' y " (E 1 - E 2 + iO)(E 1 - E 3 + i0) 



This is seen by summing terms within the columns of table (62 1 and by using the SW formula (|53|. 
The second term of Ai is a distribution with poles at E 2 — i 0, E 3 — i not pinching the i?i-axis. 
It thus vanishes to first order at E 2 = E 3 and cancels in Ai^ against the second term of Ai. We 
are thus left with 

A '' ± '° 2 "^'- £ '' U-i + i - E2 -E,± i o )' (65) 

and we conclude that Ai (+) = and Ai (-) = (2tt) 2 S(E 1 - E 2 )5(E 2 - E 3 ). The conditions E 1 = E 2 
and E 2 = E 3 are satisfied along the diagonals k\ = +k 2 resp. k 2 = +k 3 . That happens jointly and 
within the integration domain only for k\ = —k 2 = —k 3 with k\ e [fc R , & L ], whence 



as claimed. The last equality is by Eqs. (49 56) with E = k\. 

Similarly, by A n (±) (£;3 -E X ,E 2 - E 3 ) = -A i I ±) (E 2 - E 3 ,E 3 - E x ), we obtain A n (+) = and 
A„ (_) = -{2t:) 2 5{E 1 - E 2 )8{Ei - E 3 ). Hence 

An + B„ + Cb -i £ d* ^WW)" - dET(E) 2 R(E) , 



2tt 



where the last equality follows by Eqs. (49 50 1 with E = k\. 



3rd cumulant of the second kind. The cumulant is given in Eq. (40 1 as 

«(AQ) 3 » = i£d 3 i(<<^*(W3)>>+3<<^*(/iJ 2 )/ 3 >>+3<<I 1 ^* (/ 3 7 3 )» + ^*(Zi/ 2 / 3 )») • (66) 



A preliminary observation is useful. By <^> = <^*> we have (AB) = (B*A*), (T (A(ti)B(t 2 ))) = 
(JT (A(ti)* B(t 2 )*)), and likewise for higher products, T*-ordered products, and cumulants. Given 
that in the above expression the currents are self-adjoint and the ij's dummy variables, the two 
extreme and the two middle terms are so related. The missing item for establishing binomial 
statistics here is thus just 



lim \ fdM^W 3 )>>=~ \^ dET{E){l-T{E)){l-2T{E)) 



(67) 
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The computation is similar to that of the 3rd cumulant of the first kind, whence we refer to it for 
more details. By Eq. ( 35 I we have 



with 



D 



E 



f dh((i 1 ^*(T 2 i 3 ))) = D + E 

JO 



d 3 t «Ji T (J 2 / 3 )» (main term) , 



f 

Jo 

f dH((h[Q 2 ,T 2 ]))) (cont 
Jo 



(68) 



act term) . 



We apply Wick's rule (43 1 to both terms. The long-time behavior of the main term is extracted 
by Eq. (90); that of the contact term by Eq. (88 1. Expanding the commutators we so obtain 

l f* 
lim - | 



lim \ f dH ((hJ'ihh))) = Di + D„ , 

\ f dH((I 1 [Q 2 ,T 2 ]))) = E I + E lI , 
1 Jo 



lim - 

t-»+oo 



with 



Dr 



Ei 



(2tt) 3 



1 



d 3 k2iriS(E 1 - E 2 ) 



<1|J|2)<2|/|3)<3|/|1) 
E 2 - E 3 + i 



(69) 



JxK\./xR\J 



(2tt) 3 



J d 3 fc 27r5(^ - # 2 )«1|J|2X2|Q|3><3|/|1> - <1|/|2><2|7|3)<3|Q|1» . 



JxE\JxR\J 



Dji is obtained from Dj by the rule 7(23) introduced in relation with the previous cumulant. 



. (23) 

Likewise for En and Ei. In analogy with the claim (63) made there, the present one is 



Di + Ei 



1 p L 



r/'i- 



dE T(E)R(E) , D„ + E„ 



~ f dET(E)-U(E). 

J Hi 



(70) 



/ ' R 



By the same steps as those leading to Eq. (|64|, we find 
Di + Ei 



1 

(2tt) 3 



f d 3 fc ^ r I W(fc 1 ,fe 2! fc3)<l|/|2><2|7|3>W<3|7|l>, 

J c = + 



JxB\Jxl 



with the distributions 

Ti (±) = 2ni8(E 1 - E 2 



E 2 - Eo + i E 2 -E 3 + i0 



Hence this case reduces to the 3rd cumulant of the first kind ( 65 1 , which establishes the claims ( 70 1 . 
4.2 The limit of a large detector 

We will consider the situation of a detector extending over a region much larger than the Fermi 
wavelength. Clearly, the binomial distribution persists, this situation being a limiting case of the 
one dealt with before. The point though to be made is (i) that the contact terms in Eqs. ( 37|40 1 
vanish in the limit. Put differently, Matthews' time-ordered correlators reduce to ordinary time- 
ordered correlators, which alone account for the binomial distribution. Moreover, (ii) we provide 
an independent derivation of that latter fact. We shall analyze the cumulants separately, though 
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in very similar manners. The values of some integrals used along the way are collected at the end 
of the section. 

The large detector is modeled by means of scaling. Let Qo(x) be a fixed function satisfying 
( 28 1 . We choose the profile of the detector to be given by the function 

Q(x) = Q (x/l) , 

which for I ^ 1 retains that property, and consider it in the limit I — > oo. The scaling implies 

Q'( x ) = r 1 Q' (x/i), Q'(k) = cy (ik), QV(k) = r l (c%p(!k) . (71) 

2nd cumulant of the first kind, (i) We first show that the contact term vanishes in the limit 
I — » oo: 



lim - f dt x «[Qi,Ji]» 

t->00 I Jq 



I— > + 00 



0. 



The limit i — > cxd is superfluous, since is independent of t\ and in fact by (17 1 equal to 



«[<?,/]» =2 i tr(pQ'(x) 2 ) 



2i 
2^ 



dfej <1|Q' 2 |1>. 



Using (93 98) for Q' 2 instead of Q, the matrix element is seen to be a linear combination of Q' 2 (k) 
for k = 0, +2k\. They are of order Oil -1 ) by (71 1, proving the first claim. 



(ii) Let us now come to the main term: 



lim \ fd 2 t«^(f 1 f 2 )»= 7 ^ f 
t^oc t Jo (2tt) 2 J 

JxR\. 



d 2 k 



mm 



£i--E2 + i0 i^+co 2ir 



1 p 



dET(E)(l-T(E)) 



(72) 

The equality was shown in (57 59), whereas the limit is the second claim being made here. 

The regularization +i0 of the denominator only matters when Ei = E 2 , i.e. on the diagonals 
fci = ±/c2, and, once restricted to the integration domain, only for k\ = —k 2 e [fc R , k L ] (see Fig.|2|. 
Moreover, the matrix element (1|/|2) is a linear combination of Q'(+ki + fc 2 ). In the limit / —* oo 
their supports concentrate by Eq. ( [71] ) near the same diagonals, which by the same token get 
restricted to part of just one. This allows to: 

- Use the factorization E\ — E 2 + i0 = (ki — k 2 )(ki + k 2 + iO), as appropriate for k\ > 0, 
k 2 < 0. 



Select the corresponding expression (95 1 for (i\I\j) from Appendix |B] and therein neglect 
any terms vanishing on that diagonal. Hence, effectively, 

<1|/|2> = {h - k 2 )t(h : jr(k 2 )Q>(k 1 + k 2 ) . (73) 



The integrand of Eq. (72 1 so becomes 



(fci - k 2 )T{Ei)R{E 2 ) 



\Q'(ki + k 2 )\ 2 
k! + k 2 + i0 



The last factor depends on I in the way seen in Eq. (91 1 for p[x) = \Q'q(x)\ 2 = Qq(—x)Q' q (x). It 
can thus be replaced in the limit by C-5{k\ + k 2 ), where 



I 



dit 



Q'o(-y)Q' (v) 

v + iO 



(74) 
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Accepting for now that C_ = — in, we obtain the limit (72 1 by means of 2fcidfci = AE\ and by the 
earlier remark restricting k\ to [& R , fc L ]. 

3rd cumulant of the first kind, (i) We first show that the contact terms vanish in the limit 
I — > oo: 



lim - f d 2 t 0{h[Q2, J 2 ])» > , lim i f d*i [Quh]])) 

t— >00 t Jq i— H-00 t— >00 t j 



l— > + 00 



0. 



The limit £ — > oo is superfluous in the second claim, since the integrand is time-independent. It 
actually vanishes even at finite I because 

«[Q, [Q, /]]» = «[Q, [QJ]]» = 2* tr(p[Q, Q' 2 ]) = . 



The second equality is by Eq. (17 1 and the last one by the vanishing commutator. Turning to the 

1st contact term, it is convenient to use the identity ((ii[Q2, ^2])) = ^ji[Q2, ^2]))- By Wick's rule 
(42 1 and Eq. (17 1 it may then be recast as 



lim - f &H0{UQ 2 ,h}))} = f d 2 fc 

t-00 t J (2tt) 2 J 



<l|/|2><2|g' 2 |l> + <l|Q' 2 |2><2|/|l> 



E 1 - E 2 + i 

JxB\J 

By the results of Appendix |b| the numerator is a linear combination of Q' 2 (+ki + k2)Q'(+k\ + k%) 
with various sign combinations. They are of order 0(Z _1 ) by (71), proving the second claim. 



(ii) Let us now come to the main term. We showed in Eqs. (60 1 that 

' I" 



I ft . „ 

lim - d 3 t «T (7 X J 2 J 3 )» = A x + A„ , 

t— >00 r Jq 



with 



6 r d3l . <1|/|2><2|/|3><3|7|1> 

JxM\JxK\ J 



Al W J d fc (r, - r,> - / ())(£■ - zr, + / o) ' 



a = 6 f 3fe <1|/|2)<2|I|3)<3|J|1) 

11 (2tt) 3 J (^!-7; 2 + i0)(S3-^ 2 + i0) 

Jxl\JxJ 



The claim is now 
Ai 



> -1 pd^r^)^) 2 , A„ ► ~ r AET(E) 2 R(E). (75) 

Z->+00 27T l„ i^+GO 27T l„ 

17 Mr ^Mr 

It independently confirms binomial statistics, in view of TR 2 — T 2 R = T(l — T)(l — 2T). 

The computation is similar to that of the 2nd cumulant, whence we refer to the discussion 
following (72 1 for more details. We first discuss term Ai. The regularization of the denominator 
only matters when E\ = E2 or E\ = £3 and, once the integration domain is taken into account, 
only for k\ = — k 2 or k\ = — fc 3 , both along k\ e [fc a , fej. A matrix element (i ^ j) 

concentrates near the planes ki = +kj as I — > +00; and their product near the intersection: 
k\ = — k<i = — k% with fci e [A; R , fcj- This allows to: 

- Use the factorizations E\ — Ej + iO = (fci — kj){k\ + kj + iO), (j = 2, 3), as appropriate for 
k\ > 0, fc 2 , k 3 < 0. 



20 



Select the appropriate expressions for and simplify them as done for Eq. (73 1. Hence, 

we also have from (99] 97 1 



<2|J|3> = (k 2 + k 3 )(Q'(k 2 - k 3 ) - r(k 2 )r(k 3 )Q'(k 3 - k 2 )) 
<3|/|1> = (fci - k^tik^TiksjQ'i-h - kx) . 

The integrand of Ai is thus recast as 



(76) 
(77) 



(fe + *a) nEi) ® { *l + y n ^ l\ ^ (r(fe) r(fc 3 ) Q'(fc 2 - k 3 ) R(E 2 ) R(E 3 ) Q'{k 3 k 2 )) . 
(«i + k 2 + i 0)(fci + k 3 + 1 0) 

As mentioned, the expression depends on / through Q'{k) = Q' (lk); moreover, it consists of two 
terms, to each of which Eq. (92 1 is applicable by the following observation. Each term contains 
a product of distributions in the variables x = kx + k 2 , y = kx + k 3 , and x — y = k 2 — k 3 . The 
distributions correspond to 

Pi(x) = Q' Q (x) p 2 (x) = Qo(-x) p 3 (x) = Q' Q (±x) , 

where the + refers to the first, resp. second term. In the limit I — » +co, the integrand thus reduces 

to 

(ha + k 3 ) T{E X ) (C+r(k 2 ) r(k 3 ) - C^R(E 2 ) R{E 3 ))5{kx + fc 2 )<5(fc 1 + k 3 ) , (78) 

where 

y (u)Q' (-v)Q' (±(u-v)) 
(u + i0)(v + iO) 



C 4 



I 



dudv- 



Accepting for now that C + = and CL = — (27r) 2 /6, we obtain the first limit (75 1 by means of 
2k±dki = dEi and by the earlier remark restricting k\ to [fc R ,/c L ]. 

We now turn to An. In view of the integration domain the integrand of An is supported in the 
limit / — » +co near the segment k\ = —k 2 = k 3 with k\ e [fc R , fcj. We may thus: 

- Use the factorizations Ej — E 2 + iO = (kj — k 2 )(kj + k 2 + iO) (j = 1,3) as appropriate for 
k\, k 3 > and k 2 < 0. 



Select and simplify the relevant matrix elements of current as was done for Eq. (73 1. 

Hence we also have from (|97)|94|: 



<2|/|3> = (fc 3 - k 2 )t{k 3 )r{k 2 )Q'{-k 2 - k 3 ) 



<3|7|1> = (k 3 + ktXhXk^Q'ih - fci) . 
The integrand of An then reduces to 

afc + k 2 )Q'{-k 2 - k 3 )Q'(k 3 - kx) 



(79) 
(80) 



(kx + k 3 )T{E x )T{E 3 )R{E 2 )- 



(kx + k 2 +i0)(k 2 + k 3 + i0) 



In view ofQ'(fc) = Q' (lk) Eq. (92 1 may be applied with a; = kx + k 2 , y = k 2 + k 3 andx~y = kx — k 3 . 
Comparing with the derivation of the integrand (78 1 of Aj, that of An reduces in the limit to 

+ k 3 )T(E 1 )T(E 3 )R(E 2 )5(kx + k 2 )8(k 2 + k 3 ) . 



Now the second claim (75 1 follows just as the first one did from (T78 
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3rd cumulant of the second kind. By the earlier observation following ( 66 1 we only need to 
investigate 

f dh «ix^*(T 2 T 3 ))) = f dh ((hJ 1 (f 2 / 3 )» + f dh «fi[Q 2 , f 2 ]» . 

Jo Jo Jo 

(i) We first show that the contact term vanishes as I — > oo. Using ((/i[Q 2 , ^2])) = ^ji[Q2, ^2])), 
the time integral is given by Eq. (88 1. Hence, by Eq. (17), we have 



lim - \ d 2 t</i[Q 2j / 2 ])» 



2i 
2^ 



d^fe^x - £ 2 )<i|/|2><2|Qil> 



According to Eq. (95 1 and Eq. (96) with Q' 2 substituted for Q, the integrand is a linear combi- 
nation of Q' 2 (+ki ± fc 2 )Q'(+fci + fc 2 ) with various sign combinations. By (71 1 they are of order 
0(Z _1 ), proving the claim. 

(ii) We now come to the main term. We showed 



with 



Di 



lim 

t— »+co 



/2/ 3 )» = Di + D11 , 



(27T) 



2^ r fct dfci r 

J fc „ 2fcT J 



dfc. 



(2 



2i r dfci 



<1|/|-1><-1|/|3><3|7|1> 
E 1 - E 3 + i 

<1|J|-1)<-1|J|3X3|J|1) 
#3 - E x + i 



where we evaluated the delta distributions in Eq. ( 69 1 . The claim now is 

- rdET(E)R(E) 2 , D n ► ~ P d£ T(E) 2 R(E) , (81) 

/— »-+oo Z7T I,, Z^ + oo Z7T I., 



d 



which is sufficient in view of TR 2 - T 2 R = T(l - T)(l - 2T). We first consider Di. Taking the 
integration domain into account, the support of its integrand concentrates in the limit near the 
diagonal kx = — k 3 for kx e [fc R , fc L ]. This allows to: 

- Use the factorization Ex — E 3 + i = (kx — fcaXfa + k 3 + i 0) as appropriate for kx > and 
fe 3 < 0. 



Select the appropriate matrix elements of currents and simplify them as was done in Eq. ( 77 1 ; 
moreover, substituting —kx for fc 2 in Eq. (76) we obtain 



<-l|/|3> = (k 3 - kx){Q'{-kx - k 3 ) - r(kx)r(k 3 )Q'(kx + k 3 )) ; 



<1|J|-1> is given by Eq. (08). 
The integrand of Di may thus be recast as 

Q'{-kx - h) 



(k 3 -kx)T(Ex) 



(kx + k 3 + i0) 



r(kx) r(k 3 ) Q'(-kx - k 3 ) - R(Ex) R(E 3 ) Q'(k x + k 3 )) 



As mentioned, the expression depends on I through Q'(k) = Q' (lk), whence Eq. (91 1 may be 
applied with x = kx + k 3 . In the limit I — » +co the above thus reduces to 



(fc 3 - kx)S(kx + k 3 ) (C+ r(kx) r(k 3 ) - C- R(Ex) R(E 3 )) , 
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with the constants C_ and C + given by Eq. (74) resp. by 



c 4 



du 



v + iO 



Accepting on top of C_ = —iir that C + = 0, the first limit (81 I follows by 2kidk\ = dE. 

We now turn to Du. Here the support of the integrand will concentrate near k\ = fc 3 for k\ e 
[fc R , k L ] as I — » +oo. The denominator therefore factorizes as E3 — E\ + i = (fci + & 3 )(fc 3 — k\ + i 0), 
as appropriate for k\, k 3 > 0; the relevant matrix elements are given by Eqs. (48 80|) and Eq. (79 1 
with — kx substituted for &2. The integrand of Du thus becomes 



(fej + A :3 )T( J B 1 )T(S 3 )i?( J Bi) 



Q , (fci-fc 3 )Q , (fc 3 -fci) 
(k 3 -ki + i0) 



By Eq. (91 1 it reduces in the limit to C_5(ki — k 3 )(ki + k 3 )T(Ei) 2 R(Ei) , which establishes the 
second limit (81) by the aforementioned identity C_ = — iir. 



To close this section, we compute the four integrals encountered along the way. Let us generalize 
C+ to C+ = C+(0) where 

(U -v )Q'«{±{u - rV, 
v + iO 

for which we claim 

C+{u) = 0, C_(u) = -27riQ Q' (-u) . 



a 



.(u) = j dv 



Indeed, by Eq. (|5l| and Parseval's identity we have 

a 



(«) = (-*) J dv Q (-v)Q' (+(u - «)) = -2m J da; Q (x)Q' Q (+x)e+ mx . 

The first result follows by the support properties of Qo', the second is now read off and its special 
case C_(0) = — in is by $ da; Qq(x)Q' q (x) = 1/2. Let us now come to 



1 



du 



Q'o(u)C+(u) 
u + iO 



Clearly C+ = 0. By Eq. (|5l| and the SW formula (|53| we have 

I' 



du 



1 



— z 



27r»5(«))QJ,(u)C7_(u) = 2tt Jdu QoHQoQo(-") - ^-<9'o(°) 



(27T) 5 



where the last equality is by Parseval's identity followed by j da; Qq(x)Q' q (x) = 1/3, as well as by 
Q i (Q) = $dx Q' (x) = l. 

4.3 The case of a linear dispersion relation 

We shall consider the limiting case of a linear dispersion relation described in Sect. |2.2| It arises 
in the limit A — > of vanishing Fermi wavelength and thus ought to retain the binomial character 
of the transport statistics. In the model the scatterer and the detector are pointlike objects on the 
real line. We treat the cases where they are separated by I > 0, respectively coincident (1 = 0). In 
the following we briefly review the results obtained in [8] for the 3rd cumulant of the first kind, 
before extending them to the 3rd cumulant of the second kind. 
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3rd cumulant of the first kind. It was shown in [8] that the 3rd cumulant of the first kind ( 38 I 
exhibits binomial statistics, 

lim - f d 3 t«rw 2 f 3 )» = fni-r)(i-2T), 

t^+co t Jo 2n 

for both coincident and non-coincident positions of scatterer and detector. This is the counterpart 
of Eq. (34) for an energy-independent transparency T(E) = T and V = (i h — fi R . As mentioned 
in the introduction, contact terms make a non-vanishing contribution in the case I = 0. More 
precisely, we have 



lim I f d 3 i«f 7 (/ 1 7 2 f 3 )» = ^(-2T 2 )(1-T) (mam Imn! 

t-H-00 t Jq Z7T 



3 f* ^ ^ ^ 

lim - d 2 t((T(T 1 [Q 2 ,l2]))) = (1st contact term) , 

t->+00 t Jo 

If* ~ ~ V 

lim - dii «[Qi, [Qi,/i]]» = — T(l - T) (2nd contact term) . 

t->+00 £ Jq 27T 

In the second case, the main term alone contributes. In both cases a noteworthy feature is that the 
contribution of the main term remains unchanged if the time ordering is reduced as in /jT (I2I3) 
or omitted altogether. 



3rd cumulant of the second kind. As explained in connection with Eq. ( 67 I the only missing 
item in order to establish binomial statistics is 



lim \ fd 3 <«/ 1 ^(f 2 f 3 )»=f 
t-»+t» t Jq 2tt 



T(1-T)(1-2T). 



Once again, the cumulant consists of a main and a contact term, see Eqs. (68). The contribution 
of the main term is known, being insensitive to time ordering; hence we merely need to show that 
the contact term yields the complementary contribution. We distinguish between two cases: 

(i) Coincident positions, 1 = 0: From the computation of the 1st contact term in (8l we infer 



as claimed. In the second equality we used the fact that, as distributions, the odd part of (x — i 0) 2 
is ((sc - i 0)- 2 -(x + i 0)- 2 )/2 = ±m8'(x). 

(ii) Separated positions, I > 0: As mentioned in Sect. 2.2 any contact term vanishes as a 
consequence of [Qi,I{\ = 0. 



5 Comparison between different approaches 



There arc ways and means to compute cumulants of transported charge, and in particular those 
of |4j, which at first sight differ from the ones used here on various counts. The purpose of this 
section is to show that they nevertheless are fundamentally the same. 

A first difference rests on the use of Eq. pj, as opposed to Eq. (12) as above. The link is 



provided by identity (14 1, which deserves proof anyway. We recall once more that its two sides are 



to be understood as power series in A with the time ordering placed inside the multiple integrals. 
The l.h.s. equals 



e iHt e -iH{X)t 



e iAQ(t) e - 



iXQ 



W exp(i(Q(t) - Q)) = J] y -^-T{Q{t) - Q)) n ■ 



n=0 
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In order to end up with things placed as stated, we apply the fundamental theorem of calculus to 
T(Q(t) - Q)) n , rather than to (Q(t) - Q)) n : 



^m) - Q)T = f dti J- ^((Q(ti) -<?)•■• (Q(*n) - Q)) 

-i 

Jo 



t 2 =...=t„=t 



dti • ■ ■ dt„ A . . . A^((Q( tl ) _ Q) . . . (Q( tn ) _ Q)) 



We then expand the correlator in Q(ti) and — Q; the resulting second term is independent of tj and 
does not contribute to the derivative. By (13 1 this proves (14). A further proof of that identity 
also given in (8], is by comparing its two sides at each order A n . On the l.h.s. it is to be noted that 



H l (\) = -\I-i^-[Q,I] + 0(\ 3 ) 



(82) 



is not homogeneous in A; whereas the r.h.s. is expanded into ^-ordered products plus contact 
terms. It is instructive to check the case n = 2. Up to a common factor — A 2 /2 the two sides are 



r dt x dt 2 ^(/(toifta)) + r dti [Qih^ifa)] , r dhdt 2 ^*(/(ti 

Jo Jo Jo 



by (23 1 they agree. 

A further point deserving attention is as follows. For the model with quadratic dispersion the 
expansion ( |82| terminates at order A 2 and reads in first quantization 

ffi(A) = -\(pQ'{x) + Q'{x)p) + X 2 Q'(x) 2 , 



see Eqs. (16 IT I ; and with A = §dxdx'ifj(x)* A(x, x')ip(x') in second quantization 



H l {\) = -\I-i-[Q,I] 



^dx(-\j{x)Q\x) + X 2 p(x)Q'(x) 2 ^ 



(83) 
(84) 



where charge and current densities are expressed in terms of fermionic operators ip(x) as 

p(x) = $(x)*$(x) , %x) = -i${x)*${x) - ft(x)*$(x)) . (85) 



Eq. (84 1 follows from (83 1 by the commutation relation [Q, I] = [Q, I], but may also been obtained 
from H(X) as a starting point. Computations like those seen in Sect. [4] can also be performed on 
the basis of Eqs. pi [84|. However, a fact which was crucial there, but may escape notice here, 



is that the term of order A 2 in (84 1 is a commutator, as seen in ( |83| . The commutation relation, 
which states 

i[[dxp{x)Q{x)Adx'%x')Q'{x')\ = - ^ dxp(x)Q'(x) 2 , 
has a seemingly independent derivation by means of 

i[p(x),j{x')] = -6'(x'-x)p(x'), 



which in turn follows from ((85J and from [A*B, C* D] = A*{B, C*}D - C*{A*,D}B for annihi- 
lation operators A through D. 
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Appendices 



A Some limits of distributions 

In this appendix we collect some limits of distributions which are used in the main text. We 
consider distributions in one or two real variables, x and y, and parametrized by t > (or I > 0). 
In particular, we compute their limits as t — > +00. We recall the meaning of convergence for 
distributions: D t (x) —* D(x), (t — » +co) stands for 



lim dx D t (x)(j)(x) = dx D(x)4>(x) 
'^+ c0 Jr Je 



for any test function 0; say <p e <S(R), the Schwartz space 21 
Proposition A.l. We have 

p ixt _ ^ 



1 

JO 



dti e 



ixti 



1 r r 1 1 

- dtj di 2 e Mt2 = - 
* Jo Jo 

7 I dt x f dt 3 | 
t Jo Jo Jo 



— e lxt + ixt 
tx~ 2 



dt 2 e ix{tl - t2) 



t-H-oo x + iO 
i 



t-H-00 x + iO 



2tt8{x) 



(86) 
(87) 



t-> + 00 



dt-i e ixtl e iyt2 e~ i ( x+y ^ t3 



KI' d " e "") (H> 



3 e %*2 g- i(x+j/)i 3 



f^+cc (a; + i0)(ai + y + iO) 
2iriS(x) 



t— >+00 



2/ + iO 



(89) 
(90) 



Proof of (86). We have 

rt rt fOO 

lim dtte lxtl = lim d*i e ixtl e- eil = lim dh e l 

t-+aoj ^+°°Jo ^i" Jo 



The computation is justified, because the limits t — > +co and e J, interchange. This is seen by 
testing the middle expression with <f>: For either order of limits one obtains dti 0(— where 
4>(—ti) = $ R dxe zxtl 4>(x). Finally, the expression under the limit on the r.h.s. equals {(x+ie)^ 1 . □ 



As a general fact, convergence implies in the Cesaro sense; specifically lim 



Proof of 

D implies lim t _^ +00 t" 1 ^dii D tl = D. By applying this remark to (86) the claim follows 



i-»+oo Dt 



□ 



Proof of 



The integral factorizes and equals by Eq. ( 86 1 



1 e ixt - 1 e' lxt - 1 



IT 



-IT 



2 - e lxt - e 
tx* 



t— > + 00 



2tt6(x) 



The limit was established by observing that the expression on the l.h.s. is (twice) the real part of 
that in (87 1. On its r.h.s. we used the SW formula (53 1. □ 



Proof of ( 89). By the remark made earlier the average t 1 J Q dt\ may be omitted from its l.h.s. 
and the limit replaced by t\ — > +co. The remaining double integral is computed by means of 



rtl rt2 

dt 2 \ d 

Jo Jo 



dU e lyt2 e- i( - x+y)t3 



iyti 



+ 



x(x + y) y(x + y) 
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as 



X 



e i{x+y)ti _ -y 

x + y 



x(x + y) y(x + y) y 

The prefactor on the r.h.s., y~ 1 , may be given a regularization, (y + z0) _1 . It is dispensable, 
because the expression inside the brackets vanishes for y = 0, but it allows to remove them. More 
precisely, either term so obtained is a well-defined distribution; in fact a product distribution in x 
and y, respectively in x + y and y. By (86 1 the sought limit is 

1 



1 



1 



y + i0 \x + iO x + y + iO 



(x + i0)(x + y + i0) ' 



□ 



Proof of (90). The integrals on the l.h.s of (89 90 1 only differ by their domains. In fact, while in 
the first expression we demand t\ > t 2 S= £3, in the second only t% > £3 is required. This is however 
equivalent to adding the cases £2 5= t\ > £3 and t% 3= £3 > ii to the previous one. By (89) this 
yields in the limit 



1 



1 



1 



(x + i0)(x + y + iO) (y + i 0){x + y + i 0) (y + i0)(-x + iO) 



2ni S(x) 



1 



y + iO 



where we applied the SW formula (53 1 to the third term on the l.h.s. 
Proposition A. 2. For any test functions p, pi, (i = 1,2,3) we have 
p(lx) ( f p(v) 



x + i i->+oo 

Pi(lx)p 2 (ly)p 3 (l(x - y)) 
(x±i0)(y + i0) 



+ i0 



2-»+oo 



f 



dudv 



p 1 (u)p 2 (v)p3(u - V) 

(it ± iO)(v ±i0) 



S(x)S(y) . 



□ 

(91) 
(92) 



Proof of {91). We probe the convergence on any test function fa By substituting x with x/l the 
claim may be recast as 

! '% «*/()- «o» 



f 

J:-, 



da; 



0. 



x + iO v ' * ' ' " i-»+oo 
The mean value theorem yields \<fi(x) — 0(0) | < || 0' || 00 1^1 an d the sufficient bound 



dx 



p{x) 




X 


X 


10' loo 


I 



□ 



Proof of \ 92). The proof is similar to the previous one. To be shown is 

pi{x)p 2 (y)p 3 (x - y) 



J 

J 



(a; + iO)(y + iO) J-h-oo 



However, instead of using the mean value theorem, we may write 

(f>{x,y) - 0(0,0) = £01 0*0 + 2/<My) + xy<j) l2 (x,y) , 
for some further test functions fa, (i = 1,2, 12). Three terms then require estimation. The first 



one is I times 



Je 



dx p 1 {x)p{x)fa{x/l) = 0(1) , 



where p(x) = J R dy (y + zO) 1 p 2 (y)p 3 (x — y) is a bounded function; likewise for the second term. 
The third one is 0(r 2 ). □ 
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B Matrix elements of current and charge 

In this appendix we will deri ve the matrix elements of current 7 and charge Q in the quadratic 



dispersion model of Sect. |2.l| We also obtain a splitting <1|7|2> = <1|7|2)(+) + <1|7|2)^ suitable 
for Eq. d54b, i.e. 

«ibb>- <y., + <i|/|2> " 



E x -E 2 + iQ E 1 -E 2 -i0 

where (i = 1,2) denotes the Lippmann-Schwinger state and Ei = kf. 

We shall first compute (1|Q|2). Since Q(x) is supported in x > Xo, the wave-function ipk(%) of 
1 1) will matter only in that region, see Eq. ( |25[ ) . Next we will deduce (1|7|2) from Eqs. ([52j[5T}, 



from 



<1|7|2> = i{E x - S 2 )<1|Q|2>, Q(k) = {-i)- Q ' [k] 



k-iO' 



where k = +k\ + fc 2 , with signs occurring in various combinations. Finally the splitting is obtained 
by factoring k — i out of E\ — E 2 + i 0, with a sign depending on cases. 

• ki, &2 > 0. 

<1|Q|2>= T 6xfaJxjQ(x)1> ka {x)= C dxt(kJe- ik ^Q(x)t(k 2 )e ik * x 

J— oc J— oc 

= tih)t(k 2 )Q(ki - k 2 ) . (93) 

The factorization is E\ — E 2 — i = (fci + k 2 )(ki — k 2 — iO), since fci + k 2 > in the case 
considered. It is used twice: Without rcgularization, to obtain 



<1|/|2> = (fci + fc 2 )t(fciMW(fci - k 2 ) ; (94) 

and with it, to get <1|J|2>( - ) = <1|/|2> and <1|J|2>(+) = 0. In particular, by Q'(O) = 1, we 
obtain Eq. ([50]) . 

fci > 0, fc 2 < 0. 

<1|Q|2>= f dxt(kYje- tklX Q(x)(e* k2X + r(k 2 )e- lk2X ) 

J— oc 

= t(fci)Q(fci - fc 2 ) + 1Jki)r{k 2 )Q(ki + k 2 ) . 

The factorization for the second term is E\ — E 2 — i = (fci — k 2 )(ki + k 2 — i 0); but remains 
arbitrary for the first one, E\ — E 2 + iO = (fci + k 2 )(ki — k 2 — zO), since fci — k 2 here. 
Hence we have 



<1|/|2> = (fci + fc 2 )i(fci)Q'(fci - k 2 ) + (fci - k 2 )t{ki)r(k 2 )Q'{ki + k 2 ) , (95) 
and we may again set (1|/|2)(~) = (1|7|2). 

fci < 0, k 2 > 0. This case follows from the previous one by complex conjugation and by 
interchanging fci and k 2 : 



<1|Q|2> = t{k 2 )Q{h - k 2 ) + t(fci)r(fci)Q(-fci - hi) , (96) 



<1|/|2> = (fci + k 2 )t(k 2 )Q'(kx - k 2 ) + (fc 2 - fci)t(fc 2 )r(fci)Q'(-fci - fc 2 ) . (97) 
We still have (1|7|2)(~) = (1|7|2), since —i{E\ — E 2 — zO) is invariant under the move. 
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• ki, k 2 < 0. We have 

rOO 

<1|Q|2> = dx (e~ lklX + Hkije lklX )Q(x)(e tk2X + r(k 2 ) e - lk2X ) 

J— oo 

= Q(h - ha) + HM)Q{-h - k 2 ) + r(k 2 )Q(h + fca) + ^Mjr(k 2 )Q(k 2 - ki) . (98) 

The arguments of the two middle terms, + (ki + k 2 ), are non- vanishing in the case considered, 
hence their regularization by — iO irrelevant. They may thus be linked to either term 
2)(-). Using (fci + k 2 ){k\ — k 2 + i 0) = E\ — E 2 + i on the remaining two terms we have 

<1|/|2> =(fci + fc 2 )Q ; (fe 1 - k 2 ) + (fea - hWi-h - k 2 ) 

+ (h - k 2 )r(k 2 )Q'(ki + k 2 ) - (ki + k 2 )Hki)r{k 2 )Q'{k 2 - k x ) (99) 

and we may choose to set 

<1|/|2>(-) = -{k x + k 2 )Hki)r{k 2 )Q'{k 2 - k x ) . 
In particular, we find Eq. ( |56[ ). 
Acknowledgments. We thank M. Reznikov for discussions and encouragement. 
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